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Abstract 

We study the effects of non-magnetic impurity vertex correction on nu- 
clear spin-lattice relaxation rate \ jT\ of conventional s-wave superconductors 
within the Eliashberg formalism. We obtain, with a self-consistent t— matrix 
treatment of impurity scatterings, the expressions for impurity vertex function 
and nuclear spin-lattice relaxation rate. The 1/T\ is evaluated with a simple 
approximation on angular average, and found to agree in clean limit with 
the previous result that 1/(T±T) remains unrenormalized under the impurity 
vertex correction. As dirtiness is increased, on the other hand, the coherence 
peak in l/(TiT) is found to increase due to the impurity vertex correction. 
This result is discussed in connection with the experimental observations on 
conventional superconductors. 
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Nuclear magnetic resonance (NMR) relaxation rate T{ 1 measurements have been used 
extensively to study both normal and superconducting properties of various materials [|TJ. 
In simple metals, the conduction electrons provide the dominant relaxation channel for 
polarized nuclear spins, which results in the well-known Korringa law oc T. Below 
the transition temperature, T c , the superconducting gap opens and, consequently, density 
of states (DOS) increases. The former and the latter, respectively, suppresses and increases 
the relaxation rate, so that the relaxation rate is determined by the competition between 
the two: At temperatures immediately below T c , the DOS effect dominates and there is 
an increase of (TiT) -1 relative to the normal state Korringa value at T — T c . At the 
low temperature regime, on the other hand, the effect of gap opening takes over and the 
relaxation rate falls off exponentially or algebraically depending on the symmetry of the 
superconducting gap. This peak as a function of temperature is known as coherence peak, 
or Hebel-Slichter peak in NMR, whose observation was very important in establishing the 
Bardeen-Cooper-Schrieffer (BCS) theory [@-[§. 

There are several factors that affect the NMR coherence peak of {T{Ty l [|5|. The sup- 
pression of NMR coherence peak may be attributed to (a) gap anisotropy H/non s-wave 
pairing of superconducting phase |6]|7|], (b) strong coupling damping P~PX[|, (c) paramag- 
netic impurities in samples p|,|TT|], and/or (d) strong Coulomb interaction such as param- 
agnon/antiparamagnon effects |l2|Jl3ll . The non-magnetic impurity scatterings, on the other 
hand, have no influence on (T^T) -1 [|TT],|l4|,|l5j , other than the smearing of gap anisotropy. 
This is because (a) there exists a simple scaling relation between the renormalization function 
of pure (Z(u)) and impure (Z(u)) superconductors and (b) in the expression for (TiT) , 
the numerator and denominator are of the same powers in Z(u), as will be discussed in more 
detail later. Experimentally, it was observed in a series of NMR experiments on Al- [16| 



and In-based alloys [|T7| that the coherence peak in (TiT) 1 was increased as the samples 
become dirtier, or the scattering lifetime, r, is decreased. This is understood to be due to 
the increase of DOS as mentioned above as a result of gap anisotropy smearing by impurity 
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scatterings P,|18|. We point out in this paper that this may also be understood in terms of 



impurity vertex correction (IVC). We found, as will be detailed below, that in the weak cou- 
pling limit the ratio of the coherence peak to normal state Korringar value roughly doubles 
as one goes from clean to dirty limit in rough agreement with experimental observations 
P JTB|Ji7| . The enormous enhancement of NMR coherence peak reported in some organic 
superconductors, on the other hand, may be attributed to thermal fluctuations of magnetic 
flux lines pf. 

It is well known that the IVC is very important in the transport properties of metals [20 



The scattering lifetime r is changed to the transport lifetime r tr due to IVC, which appears 
in, for example, cr(0) = ne 2 r tr /m, where n, e and m are, respectively, the number density, 
charge and effective mass of the charge carriers. Because T-f" and microwave conductivity, 



lim^o have the same "plus coherence factor" [11 1, we may expect that the IVC may 
also be important for NMR relaxation rate. Maki and Fulde (MF) |21j studied this problem 
some 30 years ago assuming an infinite Fermi energy, ep. They found the impurity vertex 
function T(u,q) — > 1 in the limit of the momentum transfer q — > oo and u — > 0, and 
concluded that T± , consequently, remains unrenormalized under IVC. This conclusion, 
however, is not rigorous because Tf~ is given by an integral over the momentum transfer q, 
not by the limit q — > oo. As will be discussed later, T{~ involves integrals over k' and q of 
the product Tr G(k')G(k' + q)T(q) , where Tr and G stand for trace and Green's function, 
respectively. Since superconducting pairing occurs mainly within a narrow region around 
the Fermi surface, we can see that the main contribution to T{~ comes from the region 
\k'\ m pf and \k' + q\ ~ pp, that is, from the region Q = (q 2 + 2q ■ k')/2m ~ 0, where pp 
is the Fermi momentum. The Q ~ region and the q — > oo considered by MF should be 
distinguished. Moreover, the impurity vertex function which enters T^ 1 has a singularity as 
q — > 0, which contributes to the Q ~ region, as will be discussed. These observations call 
for a more careful analysis of the IVC effects on spin-lattice relaxation rate than the q — > oo 
of MF. In this paper, we report such a calculation. We find an explicit expression for Xf 1 
including IVC with a simple approximation about angular average which agrees with the 



approximate conclusion of MF that T^ 1 remains unrenormalized under IVC for £/£ ^> 1, 
where £ = vpr is the electron mean free path and £ = Tivp /vrA is the superconducting 
coherence length, A the gap function at zero temperature. As is reduced, the NMR 
coherence peak is found to increase due to IVC. 

The nuclear spin-lattice relaxation rate T{~ is given by |22j 

1 1 I 2 r -I 

77 = 1 - e-P» ^ \ Aq Im [X+-(u + ti,d\, (!) 

where is a form factor related with the conduction electron wavef unctions, and Q) 
is a spin-spin correlation function at frequency u and momentum transfer q, obtained 
through analytic continuation from the Matsubara function x H (iu,q). It should be cal- 
culated with non-magnetic impurities included fully self-consistently. We do this by renor- 
malizing the quasi-particle Green's function due to impurity scattering in the i-matrix ap- 
proximation and also by including the impurity vertex correction in ladder approximation, 



in accordance with the Ward identity [pjpO|. This procedure results in the Bethe-Salpeter 
equation, which is to be solved for the vertex function. The renormalized 2x2 matrix single 
particle Matsubara Green's function due to impurity scattering is given by 

Wn 2 + il + <i>l 

where Tj's are the Pauli matrices operating on the Nambu space. iW n = ip n Z n , (j) n = A n Z n , 
where p n = ir(2n + l)//3 is the Matsubara frequency, and Z n = Z(ip n ) and A n = A(ip n ) are, 
when analytically continued to real frequency, the renormalization function and gap function, 
respectively. A tilde on a variable denotes that it is renormalized by impurity scatterings, 
and £fc = k 2 /2m — /i, where \i is the chemical potential. The self-energy, S(ip n , k), due to 
impurity scatterings is then given by 

G^^ipn, k) = Gq 1 ^, k) - E(ip n , k), t(ip n , k) = nif(ip n , k), (3) 

as shown in Fig. 1(a), where n« is impurity concentration, T the t-matrix for impurity 
scattering, and Go is the Green's function for pure superconductors. We treat the non- 



magnetic impurity scattering in a self-consistent i-matrix approximation as shown in Fig. 
1(b) to obtain 
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U 



itN f 1 + U 2 



U 



T3 



w 2 + # 



iW n - <j) n Ti 



(4) 



Then there exists a simple scaling relation between pure (iW n , </>„) and impure (iW n , <p n ) 
superconductivity OJTJJ]. 



iW n = rjniW n , 4> n = r] n (p n , 
1 

Vn = 1 + 



1 ru U 2 



n, 



sin 2 5 



N ■ 



(5) 



2t^W^ + <P 2 J 2r ttN f 1 + U 2 ttN f 

where n« is impurity concentration and Sn the phase shift due to impurity scattering, U = 
7rN F V = tan^TV; where V is the impurity scattering potential energy. 

The impurity vertex function, T, is calculated in ladder approximation as shown in Fig. 
1(c). We have 

T(ip n , ip m , q) =r + J2 f(ip n )G{ip n , k')T(ip n , ip m , q)G(ip m , k' + q)f(ip m ). (6) 

%' 

This is the Bethe-Salpeter equation at arbitrary momentum transfer q and frequency 
ioJm-n — i-Pm — Wn — 27r(m — u) / (3 in ladder approximation, which we need to solve for 
T. Then T{~ can be calculated from Eq. (JJ) and 
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X+- ( iuJ ) = ^Yl p fe W + iu), 

P ip 

P(ip, ip + iu>) = ^2 Tr \p(ip + iu, k' + q)G(ip, k')T(ip, ip + iu, q) 

k',q 



(?) 



We will omit the subscripts of the Matsubara frequency for convenience. After the analytic 
continuation of iu — > uj + iS and the limit uj — > 0, Eq. (|7D is reduced to |20 



lim —Im 



X+_(cj + z5) 



de df(e) 



P(e -iS,e+ id) - Re P(e + i5,e + id) 



, 2vr de 

where /(e) = 1/(1 + e^ e ) is the Fermi distribution function. In order to solve Eq. (^j), we 



employ the following factorization as with Hirschfeld et. al. p3]. 



Then Eq. 



Lij(ip,ip + iu,q) 

Sij(ip, ip + iuj) 
Ti(ip,ip + iu,q) 

is reduced to 



- ^2 Tr TiG(ip, k')TjG(ip + icj, k' + q) 

V 



-Tr 
2 

-Tr 
2 



TiT{ip)TjT{ip + iuj) 



(9) 



(10) 



and Eq. (0) to 



(ii) 



W,1 



. .. 

The summations over fc' in Eq. (|7]) are carried out with a constant density of states, 

Sfc' = J^oo Npd^y \ /_i du, where v = cos(#), 8 is the angle between k 1 and g, and with the 

-> 2 

approximation — (k' + q) /2m— /i ~ ^r/ + <5(^), where Q{y) = (2ppqv + q 2 )/2m, which 
are standard approximations in the theory of superconductivity ||11|| . It is then straightfor- 
ward to evaluate Ly. For instance, 
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log 



x 2 + x — z5 
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^ + W„W m - y/ n y/J/y/ n y/ m , (12) 
1 



x 2 — x — i(L 



Q + £\4 + 

, x 2 + x + z<5 

±log 



x 2 — x + i5 s 



(13) 



where ^ = y'PV 2 + 2 , f n = yJW* + <%, x = q/(2p F ), and 5 S = (v"„ + f m )/e F = 
(y +y +l/r)/ep using the scaling relation of Eq. (|). ( ) represents an angular average 
over v, that is, (F) = \ j_ x dvF{y). Note that Eq. fljjD is expressed in terms of untilded 
variables taking advantage of Eq. (J5|) because its numerator and denominator are of the 
same powers in W and </>. 

Let us first consider the Born limit of 5n — ► 0. Then, T = Vt 3 , = riiTiNpV 2 , and Sij 
of Eq. (|]) becomes a diagonal matrix. In this case we can find an explicit solution T(iu, q) 
to Eq. ( |T0"D for arbitrary q and iuo. We will discuss the general case of finite 5n later. After 
somewhat laborious calculations, we obtain 
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We substitute this into Eq. ( |TT| ) to finally obtain 

K(t>m + W n W, 
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(15) 



Note that in the limit q (or x) — > 0, (/_(e — i8, e + i5)) and (/+(e — z<5, e + id>)) of Eq. Cp~3|) , 
which we need for analytic continuation as in Eq. (D, reduce to -1 + \(e F Tx) 2 + 0{x A ) and 
ieFTX 2 + C(x 4 ), respectively. The vertex functions given by Eq. flU]) , consequently, diverge 
as g — > 0, which is canceled by the phase factor of q 2 in Eq. (|TJ). We need to carry out the 
summation over the Matsubara frequency ip and the momentum q in Eq. (|15|), and then 
take analytic continuation and the limit u — > as given in Eq. (^). Neglecting IVC amounts 
to taking A = 1, for which we obtain the well known result |E3 

A 



1 ,00 9/ 

— — ; oc / rfe^— 



(ve 2 - A ; 
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Ve 2 - A 2 



°° , 9/ fe 2 + A 2 
de 



A s 



(16) 
(17) 



Eq. ( |T7D follows in the weak coupling limit. Eq. ( |i5|) may not be evaluated analytically, but is 
greatly simplified in the following approximation that treats the angular average differently: 

'/- + /- + f 2 + (W n W m + 0„0 m )/(2 v / nV / m ) 



A 



^ ' dq q 

o 



With the approximation of Eq. 
continuation of Eq. (0) 



;i +/-+/ + )a +/--/+) /' (18) 

|), we find, in the weak coupling limit, after the analytic 
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1 f» J/ f e 



oc / de— < — — + 



,2 



(4r) 2 (e 2 - A 2 



(19) 



T X T Ja de U 2 -A 2 e 2 - A 2 

We also considered the case where the impurity potential is not weak, that is, where the 
phase shift < 5n < tc/2. The resulting equation is extremely complicated. For uj — > 0, 
however, we were able to find an explicit expression corresponding to Eq. fll9|) with the 
approximation of Eq. (|18D: Eq. ( |I9|) is still valid for arbitrary phase shift. All that is needed 
is that = rijiiNpV 2 of the Born limit should be extended to 77- = -rf- , ttfi as given in 
Eq. (^). Eq. fll9"|) reduces to Eq. (|T7|) in the limit r — ► 00. This agrees with the previous 
result of MF that 1/Ti remains unrenormalized under the IVC In the limit r — > 0, on the 
other hand, 

1 r°° df fe 2 + 3A 2 1 

; de^-{——\. (20) 



T\T J a de I e 2 - A 2 

The dominant contribution to l/(TiT) comes from the region e ~ A. Consequently, 
comparing Eqs. (0) and fl20D , we can easily see that the ratio of coherence peak to normal 
state value is bigger in dirty limit than in clean limit by a factor of about 2. In Fig. 2, 
we show the normalized NMR relaxation rate, (TiT) n /(TiT) s , calculated from Eq. ([19|) , as 
a function of temperature for several values of £/£. We took the electron-phonon coupling 
constant A = 0.5 and the phonon frequency uo p h = 0.1 eV for representative values. The 
solid, dot-dot-dashed, dot-dashed, and dashed curves correspond to t/£ = 100, 10, 1, and 
0.1, respectively. The curves of£/£ = 0.1 and 100 are, respectively, almost indistinguishable 
from those of the dirty and clean limit. The NMR coherence peak is increased by a factor 
of 1.7 as expected, as is decreased from 100 to 0.1, which is in fair agreement with the 
experimental observations on Al- and In-based alloys P JT^JTTH . This observation raises an 



interesting possibility that the experimental observation of the NMR coherence peak increase 
may be understood in terms of impurity vertex correction rather than the gap anisotropy 
smearing. Experimental distinction between the two effects is highly desirable in this regard. 

To summarize, we considered the effects of impurity vertex correction on nuclear spin- 
lattice relaxation rate within the Eliashberg formalism. The relaxation rate was evaluated 



and found to agree in clean limit with the previous approximate conclusion of Maki and 
Fulde that 1/(71 T) remains unrenormalized under the impurity vertex correction. As the 
scattering lifetime is decreased, on the other hand, the coherence peak in l/(TiT) was 
found to increase due to the impurity vertex correction. The nuclear spin-lattice relaxation 
data measured on conventional superconductors as impurity scattering rates are varied may 
also be understood in terms of impurity vertex correction rather than the gap anisotropy 
smearing effects induced by impurity scatterings. 
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FIGURES 

FIG. 1. Feynmann diagrams for (a) self-energy, E, (b) t-matrix for impurity scatterings, 
T, and (c) impurity vertex function, T. The solid double and single lines, respectively, stand 
for the renormalized, G, and bare Green's functions, G . The dashed double and single lines, 
respectively, stand for the t-matrix, T, and bare impurity scatterings, Vr 3 , and the shaded 
triagle in (c) represents the impurity vertex function, T. 

FIG. 2. The normalized NMR relaxation rate, (T]T) n /(T-\T) s , as a function of T/T c . 
The solid, dot-dot-dashed, dot-dashed, and dashed lines were, respectively, calculated for 
£/£ = 100, 10, 1, and 0.1 with A = 0.5 and u ph = 0.1 eV. 
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